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Abstract 

Nanodroplets on chemically structured substrates move under the action of disjoining pressure 
induced forces. A detailed analysis of them shows that even in the absence of long-ranged lateral 
variations of the effective interface potential, already the fact, that due their small size nano- 
droplets do not sample the disjoining pressure at all distances from the substrate, can lead to 
droplet motion towards the less wettable part of the substrate, i.e., in the direction opposite to the 
one expected on the basis of macroscopic wettability considerations. 

PACS numbers: 68.08.Bc, 68.15.+e, 68.35.Ct 
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Chemically structured surfaces with their rather complex wetting properties jl| are not 
only ubiquitous in nature (e.g., beatles use them for harvesting drinking water out of fog 
[2I), they have also found important applications in ^ wide variety of processes ranging from 
ink-jet printing [3] and condensation heat transfer [4] to the so-called lab-on-a-chip concept 
for chemical analysis or biotechnology which, e.g., allows one to handle minute amounts 
of liquid containing DNA or proteins j^, B]). In view of the richness of the corresponding 
wetting phenomena and in view of their importance, they have been intensively studied both 



experimentally [3, 0, O, 1^ and theoretically 



11 



12 



tney na 

y,ii4 



15, 



16 



17 



18 



19 



20 



2l|. 



All the above mentioned applications of wettability patterns rely on the fact that the fluid 

nnnn 

placed on such a substrate will actually move towards the more wettable part [2, U, 15|, 16[ . 
However, recently we could show, that nanodroplets near chemical steps can move in the 
opposite direction [2l| . Due to the long range of dispersion forces of the van der Waals type. 



droplets residing on one side of the step perceive the different character of the other side 
even at finite distances from the step and move accordingly. Since the wetting properties 
(i.e., the equilibrium contact angles ^eq) are determined by the interplay between both the 
long- and short-ranged interactions, the force resulting only from the long-ranged part is 
not necessarily directed towards the more wettable side 21 1. 

In the present study we show that, in contrast to a droplet only near a chemical step, 
the direction of motion of a droplet positioned on top of a chemical step can depend on the 
size of the droplet and that, depending on the substrate characteristics, sufficiently small 
droplets move towards the less wettable side. For fiat droplets this effect is robust, i.e., 
independent of the detailed properties of the chemical step (in contrast to the phenomena 
involving the approach towards chemical steps [21]) and it can be stronger as compared with 
this latter approach. This is true in particular for substrates the chemical patterns of which 
are generated by modifying the short-ranged part of the liquid-substrate interaction poten- 
tial, as it is the case for many fabrication techniques (see, e.g., Ref. [10|]). We strengthen this 
prediction by mesoscopic hydrodjTiamic calculations for two-dimensional droplets (i.e., fila- 
ments) and demonstrate its experimental relevance for polystyrene nanodroplets on silicon 
oxide, a system for which the actual material parameters are available. 

We consider a droplet centered on a chemical step between two materials with coatings 
of different thickness and nature as illustrated in Fig. [TJ As in Refs. 
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22], we calculate 



the disjoining pressure (DJP) by integrating pair potentials over the substrate. In addition. 
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FIG. 1: A nanodroplet in contact with its corresponding equilibrium wetting film and located on 
top of a chemical step between two materials with coatings of different type and thickness. The 
contour and false color plot shows the vertically and laterally varying DJP for the example of 
two quarter substrates with different chemical composition but giving rise to the same equilibrium 
contact angle 0^^ = 9^^ = 90° (on the corresponding laterally homogeneous substrates) formed with 
the substrate surface. Due to the small size of the nanodroplet, its contact angles deviate from 



and d^g, respectively. The corresponding parameter values are C = 1, i? = 7, D = 1.227, 



and C = 1, = 14, = 2.251; lengths and the DJP are measured in units of b and cr/b , 
respectively; see the main text for definitions. 



in order to demonstrate the robustness of the aforementioned effect, following Refs 
Q we also use a. app.ox,„,ate DJP wKieh o. eaeh side of the ehe.ieal step equals the 
corresponding one of a laterally homogeneous substrate so that its lateral variation reduces 
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to a discontinuity at x = 0. 

A model chemical step as illustrated in Fig. [1] is composed of two different quarter spaces 
(edges) with each of them coated on its upper side with a different material. We model 
the intermolecular interactions between the liquid and the substrate as well as among the 
fluid particles via pairwise additive Lennard- Jones type intermolecular pair potentials, i.e., 
^aj^i''^) ~ ^ap^ I '^'^'^ ~ ^ap^ I ''^^ ■> whcrc M"^^^ and A*"^'^ are material parameters for the left 
(<l) and right (>) half of the substrate, respectively, and a and (3 relate to either liquid (/), 



edge (e), or coating (c) particles. Following Ref. [23|] the DJP n(x, y) of such a configuration 
can be calculated by integrating the liquid-substrate interactions over the substrate and by 
subtracting the integral of the liquid-liquid interaction potential over the volume occupied 
by the substrate. For the contribution of the left edge without surface coating occupying 
f]^ = {r G M?\x < A y < 0} one obtains 

Ii^{x,y)=AM^h2{x,y)-/\N^h{x,y), (1) 

with J„(x,y) = 4. |r - r'r'^t^V, AMe^/^ = p^Mu - pipf^'M^/'', and AiVe^/^ = pfiv, - 
piPe^^ N^/^; pi and pf^'^ are the number densities of the liquid and the edge, respectively. 
The integrals lu^x, y) and h{x, y) can be calculated analytically and can be used as suitable 
building blocks to obtain the contribution of the coating layer of thickness D'^ and of the 
right edge and its upper coating layer of thickness . With these, the DJP of the chemical 
step is given by 



n(x,2/) = n,l + n^, 

AM,^ [h^{x,y) - h2{x,y + D^)] 
-^N^ [h{x,y) - h{x,y + D^)] 
+ /^M^ /i2(x, y + D^)- /\N^h{x, y + D^) 
+ AMf [/i2(-x,y)-/i2(-x,y + D^)] 
- AATf [/6(-x, y) - h{-x, y + D^)] 

+ AMf /i2(-x, y + D"^)- AN^h{-x, y + D^), (2) 

with AMf" = pfMu - Pipf'M^/'' and AN^" = pfNu - pipf^Nf". For x ^ -00 
one has I12 7r/(45?/^) and Iq n/^Gy^) so that far from the chemical step the DJP 
reduces to 11^ and 11^ of a coated laterally /iomogeneous substrate on the left and right 



hand side, respectively. For small distances from the substrate surface, i.e., for y ^ 0, 
the leading terms (~ AM^^'^ TT/{A5y^)) in Eq. ([2]) are the terms AM^ Ii2{x,y) for x < 
and AM^ Ii2{~x,y) for x > 0, respectively. In order to have a wetting film of finite 
thickness b^^^, the disjoining pressure has to become very large at small distances from 
the substrate. Therefore, AM^^^ > is a necessary condition for the occurrence of an 
equilibrium wetting film. But the other amplitudes can be positive or negative. In order to 
simplify the presentation we only consider the case AN^^^ > as the relevant one for the 
experimental systems studied in Refs. 
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26(1 . Also by limiting ourself to coating layers 



thick compared with the wetting layer (see below) and to |AMe'^^'^| ^ AM!^^^, we neglect 



r<i/i> 



terms oc AM, 



<i/i> 
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26l |. For this system we introduce dimensionless quantities (marked 
by a star) such that lengths are measured in units of b'^ = [2AM^'^/(15 AA'^^'^I)"'^^^, which is 
the equilibrium wetting film thickness on a homogeneous flat substrate with the material 
of the left coating, and the DJP is measured in units of a/b^ where a is the liquid- vapor 
surface tension as an independent parameter. Thus the dimensionless DJP far from the edge 
(x ^ — / + oo) has the form 
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*<i/\> 

ch 



</t>- 



{y. + Dfy 



+ 



(y. + Df^)' (y. + Df^) 

with C^/^ = A^/'^b^/a, q^'^ = b^'^ /b^ (i.e., = 1), and 5^/^ = AN^'^/ANf''. The 
fully inhomogeneous and dimensionless DJP is obtained from Eq. ([2]) by replacing AMi^^'^ 



(3) 



by 45g^/^ C^/^/tt, AiV^^ by 6g^/^ C^/^/tt, and AiVe^^^ by 6g^/^ C^/^E^/^/vr. In 
order to avoid a clumsy notation in the following we drop tlie__ stars. The macroscopic 
equilibrium contact angles O^J'^ far from the step are given by 



cos<> 



Ktiy)dy. 



<i/i> 



(4) 



The wetting film thickness y^^^ minimizes the effective interface potential uj^^^{y) 
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dy' and both 6^^^ and y^^^^ can be expressed in terms of C^^^, B^^^, and 



A typical example for the DJP at a chemical step is given in Fig. [H The DJP Ilfx, y) 
obtained according to Eq. ([2]) is a continuous function of both x and y. Following Refs. [l3|, 
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FIG. 2: (a): The continuous and the sharp DJP induced lateral force /^^^ originating from the 
left and the right part of the substrate, respectively, on a droplet of height a and basal width 
2w with a = w = 15 centered on the chemical step as function of B^^^ which is varied together 
with D^^^ such that 9^^ = 9^^ = 90° (b). If the the left (right) substrate is chosen to correspond 
to I (II), / = + > 0. The force per unit length is measured in units of a. (c) and (d): 
Stokes dynamics of the same droplet on a chemical step between two substrates characterized by 
the materials parameters corresponding to points I (C^ = 1, = 1, = 7, = 1.227) and II 
(C"^ = 1, = 1, B^ = 14, = 2.251) in (a) and (b), respectively, for the continuous (c) and 
the sharp (d) DJP. The dashed line corresponds to the shape-relaxed initial state at time t = 70 
and the solid lines to t = 5000 (c) and t = 10500 (d); times are given in units of fib^/{C'^a). The 
droplets move to the right even though the equilibrium contact angles on both sides are equal. 



Iisl we also consider the discontinuous approximation 



u^hiy) for X < 
u^hiy) for X > 



(5) 



with Il'^i('^{y) defined in Eq. ([3]). In the following, we denote these two forms as "continuous" 
and "sharp", respectively. 

The DJP induced lateral force on a two-dimensional droplet (i.e., the force per unit length 
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FIG. 3: (a), (b): Values of the parameters B^^^ and D^^^ and of the contact angle 9^J^ for 
which the absolute value of the lateral force If^^'^l (originating from the left /right hand side of the 
substrate on a droplet of height a and width 2w with a = w = 5 centered on the chemical step) is 
equal to 0.705 or 0.765 (continuous) and 0.8 or 0.825 (sharp) in units of a. (c) and (d): Results 
of the Stokes dynamics for the continuous (c) and sharp (d) DJP and for the parameter sets I 
(continuous: 5^ = 6, = 1.1246; sharp: B'^ = D"^ = 1.2414) and II (continuous: B^ = 15, 
= 2.35; sharp: B^ = 15, = 2.56) in (a) and (b) for the left and the right side of the step, 



respectively (C^ = = 1 and 



1). The corresponding equilibrium contact angles with 



the substrate are indicated in the figures. The profiles correspond to times t = 55 (dashed lines) 



and 4000 (solid lines), in units of fib'^ /{C^a). From top to bottom a 



w 



15, 9.38, and 5 for 



the continuous DJP and 15, 6, and 5 for the sharp DJP. Lengths are measured in units of b"^. 
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on a liquid filament) spanning the chemical step is given by 

f = ju{x,y)n^ds = r + r, (6) 



with r as the liquid- vapor interface and the x-component of the outward surface 

normal vector on F. Changing the integration variable from the contour length ds to dy we 
obtain for the forces due to the left and right hand part of the substrate, respectively, 



/^/^ = - / n,l(a; <0)dy+ U^,{x > 0) dy, (7) 

with the droplet apex height ?/max- The first integral is the force on the left part (x < 0) of 
the droplet and the second integral the force on the right part (x > 0). For large droplets, 
i.e., for ?/max ^ Vq^^ 1 the upper limit of integration ymax can be approximated by oo. If, in 
addition, the contact lines are at a large distance (as compared to y^^^) from the step, the 
main contribution to f^^^ stems from the part of the liquid-vapor interface in the vicinity 
of the left /right contact line where T{'^?I^{x,y) ~ Tl^^{y). Therefore for macroscopic drops 
one obtains /^^^ = =F /^/^ n^^^(?/)(iy (which is independent of the droplet shape), and by 
using Eg . (HI) one recovers the well known expression / = cos 9^^ — cos 6'^ (in units of cr) 

The equilibrium contact angles Ofq^ , which determine the force on macroscopic drops, are 
given by the integrals over 'n.'^l^{y) from the equilibrium wetting film thickness to infinity 
(Eq. (j4])). The force on a nanodroplet, however, does not depend on the values of the DJP 
for y larger than its height. Therefore, one can expect that (a) the force on nanodroplets 
depends on the size of the droplets and (b) that this force can be different for substrates 
with different DJP but with nonetheless the same equilibrium contact angle. 

In order to estimate the force on a droplet spanning a chemical step — an unstable config- 
uration for which one cannot resort to equilibrium shapes — we assume a simple analytical 
expression y{x) = yo + a [1 — (x/w) J for the droplet shape, which approximates a 

parabola of basal width 2 w and height a which is smoothly connected with a wetting film 
of thickness yo. In particular for fiat droplets, the main contribution to / stems from the 
regions in the vicinity of the contact line, i.e., for relatively large distances from the chemical 
step, where 11^'''^ is well approximated by the DJP FI^''^ of a homogeneous substrate. 

For a droplet of size a = w = 15 Fig.EKa) shows /^/^ function of 5^/^ for substrates 
with C^/^ = 1 and D^/^ chosen such that 9^^ = 9^^ = 90°. Evidently, the droplet is too 
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small to justify the sharp approximation, but in both models, is monotonically decreasing. 
Although according to the above parameterization y{x) the contact angles O^J'^ = 90° are 
fixed, the force depends on the substrate parameters such that due to + /'^ 7^ a droplet 
on top of a step between two sides with equal 6eq but different chemical composition is driven 
to one or the other side of the step. A macroscopic drop does not move at all in this case. 
The chosen parameter values (C^ = and = = 1) correspond to two substrates 
with equal Hamaker constants. Therefore, to leading order in the distance from the step, 
the lon g-ra nged lateral force on a droplet sitting next to (and not on top of) a chemical step 



is zero [2l|]. A careful analysis of the force on a slightly displaced droplet indicates, that 
the droplet stops when the trailing edge reaches the contact line. At this point the force 
changes sign. 

This behavior is confirmed by mesoscopic hydrodynamic calculations based on the 2D 
Stokes equation for an incompressible liquid exposed to either the continuous or the sharp 
DJP (see Figs. [2](c) and (d)). This dynamics also takes into account the shape relax- 
ation of the droplet. In dimensionless form the continuity and Stokes equations read 
V ■ u = and V^u = Vp, respectively, where u = {ux,Uy) is the velocity vector 



and p is the pressure 21j. With the viscosity /i, we use as the velocity and the time scale 
a / jj, and fib^/{C^a), respectively. Lengths and pressure are expressed in units of and 
<j/b^, respectively. We have solved theses equations numerically with a standard biharmonic 
boundary integral method [28]. A no-slip boundary condition has been employed for the 
impermeable liquid-solid interface and a no-fiux boundary condition was imposed at the 
left and right end of the system. At the liquid-vapor interface the tangential stresses are 
assumed to be zero and the normal stresses are balanced by surface tension and the DJP 



28|, 



29| leading to 

n ■ r ■ n = -K + n, (8) 

with the local mean curvature k, the stress tensor r, and the surface normal vector n 
pointing out of the liquid. Here, the parameterization of the droplet shape used in the force 
calculation was used as the initial configuration. 

In Fig. [2] we show that the net force on a droplet can be non-zero even for equal macro- 
scopic equilibrium contact angles on both sides. Since the contact angle as well as the force 
vary monotonically with the parameters, the net force is not necessarily directed towards 
the more wettable side of the step. As an example Figs, and (b) show various values 



of and D^/^ for which the absolute value of the lateral force on a nanodroplet (with 

dimensions a = w = h) due to left/right side of the substrate is equal to |/"^^^| = 0.705 
or l/^/^l = 0.765 (continuous) and |/^/^| = 0.8 or \f^'^\ = 0.825 (sharp). As function 
of these products 5^/^ and D'^^^, Ofq^ changes accordingly. For both values of f^^^ one 
can find system parameters such as the ones denoted as I and II so that the force from 
the less wettable side I is larger than from the more wettable side II. Consequently, for 
such systems a droplet moves in the direction opposite to the one expected on the basis of 
the difference of contact angles. As in the example shown in Fig. [21 the force analysis for 
displaced droplets indicates that, independent of the droplet size, the droplets will stop their 
motion once the trailing contact line gets close to the chemical step. Both behaviors are also 
confirmed by the Stokes dynamics shown in Figs.[3](c) and (d): a small droplet of dimensions 
a = w = b moves towards the less wettable side, while a droplet of size a = w = 15 already 
behaves like a macroscopic droplet in the sense that it moves towards the more wettable 
side. For the system parameters chosen for Fig. [3] the critical size, for which the direction 
of motion in the Stokes dynamics changes sign, is found to be a = w = 9.5 and a = w = Q 
in the case of the continuous and the sharp DJP, respectively. 

Although the phenomenon described in Fig. [3l i.e., the motion towards larger contact 
angles, is not generic, it can be expected to appear in the extensively studied experimental 
system polystyrene (PS: a = 30.8 mN/m, fi = 1200 Pas 24]) on a silicon waver covered by 
SiO ^25]. The experimental study reported in the Ref. 25] indicates that the DJP of this 
system, which is shown in the inset of Fig. IH^a) for two different thicknesses of oxide layers 
SiO (2.4 nm and 191 nm), can be represented by an expression of the form given in Eq. ([3]). 
The corresponding equilibrium contact angles are 6^^ = 7.2° (C^ = 0.054, B"^ = —5.91) and 
6*^ = 7.6° {C^ = 0.0155, = —5.91). A macroscopic droplet sitting on the chemical step 
should therefore move towards the side with the thinner oxide layer. The force calculation 
for droplets of varying sizes shows, however, a change of sign at a droplet height of about 
7.5 nm for both the sharp and the continous DJP. In agreement with the force calculation, 
the Stokes dynamics calculations predict that a droplet of height a =5 nm moves towards the 
less wettable substrate. In the Stokes dynamics the droplet stops once its trailing contact 
line reaches the step. It covers about 80 nm in roughly 1100 s (6^ = 0.84 nm), which should 
be experimentally observable. 

In summary our study of the dynamics of droplets spanning chemical steps has revealed 
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FIG. 4: (a) The DJP induced force of per length on a PS droplet centered on a chemical step 
composed of a Si wafer coated with oxide layers of two different thicknesses (I: = 2.4 nm, and 
II: = 191 nm) as a function of the droplet height a. The ratio between the height a and the 
width w in the ansatz for the droplet shape is chosen to be tan(a) with a = 7.5°, i.e., in between 
the equilibrium contact angle 7.2° and 7.6° on the left and the right hand side, respectively. For 
a ~ 7.5 nm the force changes sign. The inset shows the DJP on both sides, (b) Stokes dynamics 
of a PS droplet with an initial height a =5 nm on a chemical step composed of coatings of type I 
on the left side and of type II on the right side. In (b) the profiles correspond to times i = 0.1 s 
(dashed lines) and to t = 1050 s and 1200 s (solid lines) for the continuous (top) and the sharp 
(bottom) DJP, respectively. 



a qualitatively different behavior of nanodroplets as compared to macroscopic droplets 
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12l |. While the dynamic behavior of macroscopic droplets is solely determined by 



the difference of the equilibrium contact angles, the direction of motion of nanodroplets 
depends on their size and on details of the composition of the substrate. Nanodroplets can 
move towards one side of a chemical step under conditions for which macroscopic droplets 
either do not move at all or move to the opposite side of the step. Since most actual 
substrates exhibit nanoscale chemical heterogeneities, our results provide a basis for a better 
understanding of the dynamics of wetting phenomena, in addition to possible applications 
for handling, controlling, and guiding liquids in emerging nanofluidic devices. For instance, 
because the size, for which the direction of motion changes sign, also depends on the chemical 
composition of the droplets, this size selectivity has potential applications for droplet sorting 
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by size or by composition in nanoscale fluidic devices. 
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